Phantom energy accretion and primordial black holes evolution in Brans-Dicke theory 
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In this work, we study the evolution of primordial black holes within the context of Brans-Dicke 
theory by considering the presence of a dark energy component with a super-negative equation of 
state called phantom energy as a background. Besides Hawking evaporation, here we consider two 
type of accretions - radiation accretion and phantom energy accretion. We found that radiation 
accretion increases the lifetime of primordial black holes whereas phantom accretion decreases the 
lifespan of primordial black holes. Investigating the competition between the radiation accretion 
and phantom accretion, we got that there is an instant during the matter-dominated era beyond 
which phantom accretion dominates radiation accretion. So the primordial black holes which are 
formed in the later part of radiation dominated era and in matter dominated era are evaporated at 
a quicker rate than the Hawking evaporation. But for presently evaporating primordial black holes, 
| radiation accretion and Hawking evaporation terms are dominant over phantom accretion term and 

, hence presently evaporating primordial black holes are not much affected by phantom accretion. 
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I. INTRODUCTION 

Einstein's formulation of General Theory of Relativity(GTR) [l[ in 1916 takes gravitational constant(G) to be a 
' time-independent quantity. It is a pure tensor theory of gravity. Following Einstein's lead, many scalar tensor theories 
have been developed as the extensions of GTR. In all these theories G is a time dependent quantity. Among them 
^Q ' Brans-Dicke(BD) theory Q is the simplest one. In BD theory the gravitational constant is set by the inverse of a 
time-dependent scalar field which couples to gravity with a coupling parameter lo. GTR can be recovered from BD 
theory in the limit cj — > oo 3]. BD theory also admits simple expanding solutions [H for scalar field (f>(t) and scale 
factor a(t) which are compatible with solar system observations (5|. BD theory is also sucessful in explaining many 
cosmological phenomena such as inflation Q, early and late time behaviour of the Universe Q, cosmic acceleration 
and structure formation cosmic acceleration, coincidence problem 0, [13] and problems relating to black holes 

Primordial Black Holes (PBHs) could be formed in the early universe due to various mechanisms, such as inflation 
[li], [IB] i initial inhomogeneities [T(|, phase transition and critical phenomena in gravitational collapse [l7i [l8j. bubble 
collision [T|| or the decay of cosmic loops [2(| The formation masses of PBHs could be small enough for them to 
have evaporated completely by the present epoch due to Hawking evaporation |2l| . Early evaporating PBHs could 
account for baryogenesis [22|,[23[ in the universe. On the other hand, longer lived PBHs could act as seeds for structure 
formation [24j and could also form a significant component of dark matter [25j j . 

The finding of SN la observations [26( that the Universe is currently undergoing accelerated expansion constitutes 
the most intriguing discovery in observational cosmology of recent years. As a possible theoretical explanation, it is 
considered that the vacuum energy with negative pressure having equation of state v — IP termed as dark energy, is 
responsible for this acceleration. Some works [13, HH have raised the possibility that the equation of state parameter 
7 may be less than —1, which is known as phantom energy in literature. A peculiar property of cosmological models 
with phantom energy is the possibility of a Big Rip [2t| : an infinite increase of the scale factor of the Universe in a 
finite time. In the Big Rip scenario, the cosmological phantom energy density tends to infinity and all bound objects 
in the Universe are finally torn apart up to the subnuclear scales. The Present data on distant supernovas [3fJ showed 
that the presence of phantom energy with —1.2 < 7 < —1 in the Universe is highly likely. It has been already 
acknowledged that, being such an exotic physical species, the phantom energy may change the accretion regime of 
black holes. 

In this work, we study evolution of PBH within a general phantom energy scenario. Taking phantom energy 
accretion along with Hawking evaporation and accretion of radiation, we show the PBH evolution in different cosmic 
era and discuss about the presently evaporating PBH. 
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II. PBHS IN BRANS DICKE THEORY 



For a spatially flat(fc = 0) FRW universe with scale factor a , the Einstein equations and the equation of motion 
for the JBD field $ take the form 
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The energy conservation equation is 



P- 



|(£)(l + 7 )p = (4) 



on assuming that the universe is filled with perfect fluid describrd by equation of state p = 7/?. The parameter 7 is 
■| for radiation dominated era(i < t e ) and is for matter dominated era(i > t e ), where time t e marks the end of the 
radiation dominated era w 10 11 sec. 
Now equation Q gives 

fa- 4 (t < t e ) 
P0C \a- 3 (t>t e ) 

Barrow and Carr [3l[ have obtained the following solutions for a and G for different eras, as 

Go(^) (t>t e ) 

where to ~ is the present time, Go ~ is the present value of G, and n is a parameter related to u, i.e., n = 



and 



Since solar system observations 32j require that to be large (oj > 10 4 ), n is very small (n < 0.00007) 



From our previous work ll| , we know that if we consider Hawking evaporation and accretion of radiation symul- 



taneously, then the rate at which primordial black hole mass changes is given by 

AW --2^g4^ + 16 - G2M2/ ^ (7) 
where an is the Stefan-Boltzmann constanat and / is the accretion efficiency. 

In a Universe also filled with phantom energy the accretion of such exotic component should also be taken in 
account. Babichev, Dokuchaev and Eroshenko [27| have worked out a differential equation for a black hole accreting 
phantom energy only and found that phantom energy accretion decreases the overall black hole mass. The rate of 
accretion of phantom energy is given by 

M ph = 16irG 2 M 2 [p ph +p{p ph )} (8) 

But for phantom energy, p(p) = 7/? with 7 < — 1. So above equation takes the form 

M ph = 16^G 2 Af 2 (l + 7 )p ph (9) 

Considering the radiation accretion and Hawking evaporation term along with new phantom energy accretion term, 
one can get the PBH evolution equation as 

Mpbh = - 25 a 6 ^ 3 Q l M2 + l6TrG 2 M 2 fp r + l6nG 2 M 2 (l + j)p ph (10) 

This equation is not exactly solvable but we solve it by using numerical methods for different cosmic era. 
In our calculation, we have used 7 = —1.1. 
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III. EVOLUTION OF PBH IN RADIATION-DOMINATED ERA 
A. Accretion of radiation 

The equation for accretion of radiation is 

Mrad = 16nG 2 M 2 fp r (11) 

lergy density whii 

equation can be written as 



where p r is the radiation energy density which varies with scale factor as p r = — ^) . Using p® — f^Pe, above 



M rad = mnG 2 M 2 f Pc n° r (— ) 4 (12) 
Va / 



Again using equations §5$ and ©, we get 

/„ x 2 " > ■ / ' \ 2 /t e \<-- s+i " > " 

By inserting numerical values of different quantities i.e. Gq = 6.67 x 10~ 8 dyne-cm 2 /gm 2 , p c = 1.1 x 10 -29 gm/cm 3 , 



M rad = 16^ 2 (^)^M 2 /^(l)^(^)-™ /J (13) 
of different q 

t = 4.42 x 10 17 sec along with w 10" 5 , one can find 



/tn\ n M 2 

M rad =l.MfG Q {J) — (14) 

On integration above equation gives 

j#- w 

1-1.94/(^-1) 

where Mi is the initial mass of PBH formed at time U. 
For large time t, the above equation asymptotes to 

M- 

M = — ± (16) 

1 - 1.94/ V ' 

Thus for accretion to be effective, 

0.515 (17) 

The radiation accretion of a particular PBH having initial mass Mi = 10 10 gm is shown in figure-1, which indicates 
mass of the PBH increases with accretion efficiency. 

For simpicity, all graphs are plotted in logarithmic scale. 

B. Accretion of phantom energy 

Accretion due to phantom energy is govern by the equation 

M^ = 16^G 2 M 2 (l + 7 )p p/l (18) 
But phantom energy density varies with scale factor as 



Pph = h f | ( — ) (19) 
1 +7 Van/ 



So phantom energy accretion equation becomes 

::i i + - ] 

Again using equations ([5]) and dU), we get 



M ph = -16irG 2 M 2 p° p J— ) 3(1+7) (20) 
van/ 



Mph = -167 



/f„\2n / / N -3(l+7)/2 // \ -(2-n)(l+7) 

*<d M n Mh) © (21) 

By solving equation (|2T|) . Phantom energy accretion of a particular PBH having Mj = 10 10 gm is shown in figure-2. 
This figure indicates phantom energy accretion is ineffective in radiation dominated era. 
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FIG. 1: Variation of PBH mass for diffent accretion efficiencies as / = 0.1, 0.2, 0.3 

phantom energy accretion 
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FIG. 2: Variation of PBH mass due to phantom energy accretion 



C. Complete evolution equation 

In radiation dominated era, the complete rate of change of PBH mass is given by 

MPBH 



nr, 1 /t\-2/t \ ((-8+4ra)/3)-2ra 

t \-3(l+7)/2/t p \ -(2-n)(l+ 7 )-2n 



-167rGgM 2 /o c n! 



(22) 



The complete evolution of a particular PBH having Mj = 10 10 <7m is shown in the figure-3. It is clear from the 
figure that evaporation time of a PBH delays due to accretion of radiation. More is the accretion, more delay is the 
evaporation. 
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FIG. 3: Variation of evaporating time of a PBH in radition dominated era for / = 0, 0.2, 0.3 



IV. EVOLUTION OF PBH IN MATTER-DOMINATED ERA 



In this era, radiation accretion equation (|13|) takes the form 

/ t \ ((-8+4n)/3)-2n 

M rad = WnG 2 M 2 fp c n° r (-) 

Integrating above equation, we get 
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M(t) = M; 

For large time t, this equation asymptotes to 

M(t) = Mi 
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In matter dominated era, phantom energy accretion equation (|2ip takes the form 

t X -(2-n)(l+ 7 ) 



M, 
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Integating above equation, we get 
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1 



16 



(2-n)|l+7|+2-n 



(2 - n)|l + 7| +1 -2n 
For large time t, this equation gives 
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(26) 



(27) 



(28) 



The variation of PBH mass with time due to phantom energy is shown in figure-4 which indicates mass of the PBH 
decreases due to phantom energy accretion. 
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phantom energy accretion 



1.0 x 10 



7.0x10 



bO 48 

c 5.0x10 



3.0x10 



10' 



10' 



10 17 10 
t in sec 




FIG. 4: Accretion of Phantom energy for ti = 10 sec is shown in the figure 
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FIG. 5: Phantom accretion (Red) and Radiation Accretion (Blue) curves for ti = 10 11 sec are shown in the figure 



The variation of PBH mass for phantom energy accretion and for radiation accretion is shown in figure-5. From 
this figure, it is clear that radiation accretion increases the PBH mass whereas phantom energy accretion decreases 
the PBH mass. 

Crossover from radiation accretion to phantom accretion occurs at a time t eq which can be calculated by equating 
the magnitude of two type of accreting masses. So from equations (|25|) and (|28p. one can get 



eg 



ti x 



3{(2-n)|l+ 7 | + l-2n} 



5 + 2n 



(2-n)|l + T| + l-2, 



y/fi? 

ph 



(2-^)11+71 + 1-27, 



(29) 



t eq is the time at which cross over from radiation to phantom occurs. This crossover time t eq increases with increase 
in accretion efficiency which is shown in figure-6. 

In comparison with matter dominated era PBHs, rate of radiation accretion is much larger and rate of phantom 
accretion is much smaller for radiation dominated era PBHs. So for these PBHs the crossover time from radiation to 
phantom (t eq ) comes much later and hence some radiation dominated era PBHs are completely evaporated without 
feeling the phantom energy accretion. 
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FIG. 6: Crossover from radiation to phantom accretion for t, = 10 sec is shown for different accretion efficiencies 



In matter dominated era, the complete rate of change of PBH mass is given by 
a H 1 



M 



PBH 



256tt 3 G 2 M 2 
Solving equation (1501) . we construct the Table- 1 



+ 16 , G s(i) 3 "Mv[/n;(i) 



t \ (-8+4n)/3 
to ■ 



ph 





(2-„)(l+ 7 )- 







(30) 



U = 10 11 sec 


/ 


ievap 


(j^evap^ph 





3.333 x 10 118 sec 


4.445 x 10 43 sec 


0.2 


7.379 x 10 118 sec 


4.445 x 10 43 sec 


0.4 


2.181 x 10 119 scc 


4.445 x 10 43 sec 


0.6 


1.211 x 10 120 sec 


4.445 x 10 43 sec 


0.8 


1.007 x 10 121 sec 


4.445 x 10 43 sec 



TABLE I: The evaporating times of the PBHs which are created at t — 10 11 sec are displayed for several accretion efficiencies 
for both cases without phantom accretion (t evap ) and with phantom accretion (t evap ) p h- 

From Table- 1, we found that due to phantom energy accretion PBHs evaporated at a quicker rate than Hawking 
evaporation and their evaporating time is independent of radiation accretion. 



EVOLUTION DYNAMICS FOR PRESENTLY EVAPORATING PBH 



In this section, we discuss about the PBHs whose evaporating time is to- Solving equations ([22)) and (|3H)) numerically, 
we construct the Table-2 for presently evaporating PBHs. 

It is clear from the Table-2 that phantom energy accretion does not affect lifetimes of presently evaporating 
PBHs because crossover time t eq from radiation to phantom is much larger than evaporating time. For same rea- 
son, the PBHs which are completely evaporated by present time are not influenced by the presence of phantom energy. 

Now we calculate the 7-ray constraint on PBHs by assuming observed 7-ray background arises due to presently 
evaporating PBHs. 

The fraction of the Universes' mass going into PBHs at time t is given by [l6[ 



PBH 

Th7 



(*)- 



(31) 
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tevap = to = 4.42 x 10 17 sec 


/ 


Mi 


(Mi) ph 





2.3669 x 10 15 gm 


2.3669 x 10 15 gm 


0.1 


1.908 x 10 15 gm 


1.908 x 10 15 gm 


0.2 


1.449 x 10 15 gm 


1.449 x 10 15 gm 


0.3 


0.989 x 10 15 gm 


0.989 x 10 15 gm 


0.4 


0.530 x 10 15 gm 


0.530 x 10 15 gm 


0.5 


0.714 x 10 14 gm 


0.714 x 10 14 gm 



TABLE II: The formation masses of the PBHs which are evaporating now are displayed for several accretion efficiencies for 
both cases without phantom accretion (Mi) and with phantom accretion (Mi) p h- 

where £lpBH{t) is the density parameter associated with PBHs formed at time t, z is the redshift associated with 
time t, Qr is the microwave background density. 

If M„ be the mass of the presently evaporating PBH, then initial mass fraction for that PBH becomes [ll[ 

The variation of /3(M») with / drawn from variation of M* with / is shown in the Table-3. The bound on /3(M*) is 
strengthened as / approaches its maximum value. But the bound is independent of phantom energy. Case is same 
for all other observed astrophysical constraints which arises due to completely evaporating PBHs. 
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2.3669 x 10 15 gm 


5.71 x 10" 26 


0.1 


1.908 x 10 15 gm 


5.13 x 10~ 26 


0.2 


1.449 x 10 15 gm 


4.46 x 10" 26 


0.3 


0.989 x 10 15 gm 


3.69 x 10~ 26 


0.4 


0.530 x 10 15 gm 


2.70 x 10" 26 


0.5 


0.714 x 10 14 gm 


0.99 x 10" 26 



TABLE III: Upper bounds on the initial mass fraction of PBHs that are evaporating today for various accretion efficiencies /. 



VI. CONCLUSION 

In this work, we study the evolution of primordial black holes within Brans-Dicke theory by considering the presence 
of a dark energy component with a super-negative equation of state called phantom energy as a background. Along 
with Hawking evapoartion, here we consider two type of accretions - radiation accretion and phantom energy accretion. 
We discuss each accretion term separately and then add all accretion and evaporation terms to study the complete 
evolution of PBH for different cosmic era. We found that in radiation dominted era phantom accretion is ineffective. 
But in matter dominated era, radition accretion increases the lifetime of PBHs whereas phantom energy accretion 
decreases the lifespan of PBHs. Investigating the competition between the radiation accretion and phantom accretion, 
we got that there is an instant during the matter-dominated era beyond which phantom accretion dominates over 
radiation accretion. So the PBHs which live beyond this transition time are affected by phantom energy accretion 
and evaporated at a quicker rate than their Hawking evaporation. Mainly those PBHs which are formed in the later 
part of the radiation dominated era and in matter dominated era are influenced by phantom energy. But for the 
presently evaporating PBHs, the radiation accretion and Hawking evaporation terms are dominant over phantom 
accretion term. For these PBHs crossover time from radiation to phantom comes much later than their evaporating 
time. Hence presently evaporating primordial black holes are not affected much by phantom accretion. Case is same 
for all primordial black holes which are evaporated till now. Thus, all observed astrophysical constraints on primordial 
black holes remain unaltered in the presence of phantom energy. 
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